To Appear in the Astrophysical Journal 

Preprint typeset using I4TeX style emulateapj v. 08/22/09 



o 
o 

(N 

Oh: 

00 



< 

6 



> 

(N 

0^ 
O 

o\ 

O 



ON THE TIMESCALE FOR STAR FORMATION IN GALAXIES 

Tony Wong 

Astronomy Depailment, University of Illinois, Urbana, IL 61801 
Received 6 Aug 2009; Accepted 8 Sep 2009 

ABSTRACT 

The timescale for star formation, a measure of how quickly neutral gas is being converted to stars, is con- 
siderably longer than typical dynamical timescales associated with a galactic disk. For purposes of modeling 
galaxy evolution, however, it would be extremely attractive if the star formation timescale were proportional 
to an easily derived dynamica l timescale. W e compare estimates of the star formation timescale within nearby 
galaxies, based on the work of' Lerov et al.l (1200 8i) and existing BIMA SONG CO data, with three simple forms 
of the dynamical time: the orbital time, the free-fall time at the midplane density, and the disk Jeans time (the 
growth time for gravitational instabilities in a disk). When taking into account the gravity of the stellar disk in 
an approximate way, all three timescales show correlations with the star formation timescale, though none of 
the correlations can be accurately described as linear. Systematic errors in estimating appropriate gas masses 
and the stellar velocity dispersion may obscure an underlying correlation, but we focus instead on a model 
where the timescale for H2 formation from H I is decoupled from the timescale for star formation from H2. The 
Jeans time correlates well with the first of these timescales, but the relationship is still non-linear, and requires 
a characteristic GMC lifetime that increases toward galaxy centers. 
Subject headings: galaxies: ISM — galaxies: evolution — stars: formation 



1. INTRODUCTION 

In recent years increasing attention has been focused on the 
question of how the star formation rate (SFR) in disk galaxies 
is determined, as this remains one of the key parameters dif- 
ficult to specify ab initio in cosmological simulations. By far 
the most common prescription is the "Schmidt law", which 
posits that the SFR per unit area (or volume) scales with the 
gas surface (or volume) density raised to a fixed power, e.g. 



■^SFR ■ 



(1) 



as suggested by lSchmidl(ll959h and first demons trated empir- 
ically for external galaxies by iKennicuti (Il989h . (Hereafter 
we use to refer to the power law relating surface densi- 
ties and n for that relating volume densities.) Since the sur- 
face density S is simply the projection of volume density p 
through the disk, given a constant scale height for gas and 
star formation one should obtain the same power law indices 
for surface and volume densities. This has led to the sug- 
gestion that the observed index of ( Kennicut3ll998l) 
can be understood if the SFR per unit gas mass, p*/pg (here- 
after the star formation efficiency or SEE), is proportional 
to the free-fall rate, ojff (x (GppY^^. As is well known (e.g., 
IZuckerman & EvansI Il974^ . the actual proportionality con- 
stant must be much smaller than 1 in order to account for 
the observed inefficiency of star formation (i.e., the relatively 
long timescale for consuming aU interstellar gas compared to 
the free-fall time). 

On the other hand, the assumption of a constant gas scale 
height is unlikely to be correct; the H I disks of galaxies are 
observed to flare outwards at large radius (e.g.. Oiling 1996), 
while the vertical gas ve locity dis persion appears to be 
fairly constant (Shostak & van der K ruit 1984). For a purely 
gaseous disk in gravitational equilibrium, the scale height can 
be approximated as hg = c^/ (ttGE^) (Ivan der Kruit & Searld 
11981b . and the vertical crossing time hg/cg=Cg/(TTGT,g) is 
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analogous to the free-fall time. If the SEE indeed scales with 
the free-fall rate, one therefore expects a power-law index 
N=2 for surface density. 

The widely adopted N=l.4 Schmidt law of IKennicuti 
(1998) is based on global disk averages augmented with spa- 
tially resolved images of infrared-luminous starbursts. In gen- 
eral, observational tests of the Schmidt law within galaxies 
have been limited by the scarcity of spatially resolved CO 
and H I maps of nearby galaxies. Wong & Blitz (12002b in- 
vestigated the azimuthally averaged star formation law in a 
sample of seven CO-luminous galaxies from the BIMA Sur- 
vey of Nearby Galaxies (BIMA SONG). Using their favored 
model for correcting Ha emission for extinction, they esti- 
mated N 1 .75 ± 0.2 5 when compa r ing the surface densities 
of total gas and SFR. iBoissier et al.l (1200 3*) employed some- 
what different assumptions for inferring gas and SFR surface 
densities, but arrived at a similar result of » 2 for asam- 
ple of 16 galaxies. On the other hand, iHever et aTl (120041) 
found w 3.3 in M33, a consequence of the relatively flat ra- 
dial gas profile in this H I-dom inated galaxy. More recently, 
IHgiel et al. (200i) and lLeroy e t al. (2008) have exploited im- 
provements in H I, CO, and infrared imaging to examine the 
Schmidt law across a wider range of galactocent ric radii and 
in bot h CO-luminous and CO-faint galaxies. iBigiel et aTl 
(I2OO8I) find an average index of N 1 .85 ± 0.70 for a sam- 
ple of seven gal axies, with a facto r of ~2 scatter about the 
mean relation. Lero y et aP (120081) identify two regimes, an 
H2-dominated regime (in the inner disks of massive spirals) 
where the SEE is roughly constant, and an H I-dominated 
regime (in outer disks and dwarfs) where the SEE shows a 
strong positive correlation with the stellar surface density, S,. 

The high degree of variability in the observed power law 
index A^, and the lack of any clear correlation between the 
SFR and H I surface density in spatially resolved studies 
(see also iKennicutt et al.| 120071) . have led many authors to 
conclude that the Schmidt law is largely a property of the 
molecular gas, since restricting the comparison to CO fluxes 
yields a narrower range of values for A^ ranging from 1-1.4 



(iWong & BUtzl I200I iHever et all l2004t iBigiel et all llOOSh . 
Indeed the densest molecular gas, traced by HCN 7=1-0 emis- 
sion, appears to be converted into stars on a roughly constant 
timescale of ^10* ^ yr (Gao & Solomon 2004; Wu et al. 20051). 
iKrumholz & Thompson! (12007 ) argue that this timescale is 
proportional to the free-fall time at the critical density of the 
HCN line, which is higher than the median density of the ISM: 
thus all HCN-emitting regions tend to have this density and 
free-fall time. However, focusing solely on high-density gas 
leaves unanswered the question of how the H2 or "dense" gas 
content is determined, and whether galactic-scale processes 
(involving one or more dynamical timescales) are involved. 
Furthermore, if the gas traced by CO emission is part of a 
continuous range in density, one might reasonably expect that 
the star formation rate can be predicted from the average gas 
density, including the atomic gas. 

Our focus in this paper is on obtaining simple estimates of 
three dynamical timescales for a galactic disk composed of 
gas and stars, and comparing these estimates with recent data. 
Although none of these timescales appear to yield a satisfac- 
tory prediction for the SFR, the inclusion of the stellar com- 
ponent appears to be important in reproducing some of the 
observed trends. Our estimates are based on azimuthal aver- 
ages, and thus may fail to reflect the appropriate local dynam- 
ical timescales. Moreover, the difficulty of measuring stellar 
and gaseous velocity dispersions, and the frequent presence of 
a stellar bulge in regions of bright CO emission, will compli- 
cate the evaluation of the free-fall time. Thus, this approach 
must be tested and refined in nearby systems where star for- 
mation and GMCs can be spatially resolved. 

2. THREE DYNAMICAL RATES 

The fundamental dynamical rates for a galactic disk are 
the orbital rate Vl and the Jeans rate lo^ = ttGS/c,, where S 
is the mass surface density and Cj is the sound speed (typ- 
ically replaced by an "effective" sound speed that includes 
the contributions of turbulent and magnetic pressure). The 
two are related via the axisymmetric instability parameter 
Q = Kc,/(^GI]): 

K 17 
= t: oc 1 

Q Q 

where we have assumed k oc 17 (e.g., a flat rotation curve). 
The Jeans rate is the natural growth rate for gravitational in- 
stability in a disk, although it should be noted that axisymmet- 
ric (ring) instabilities will be suppressed if Q> 1, and even for 
( 2 < 1 their growt h rate is suppressed by a factor of \/\ — Q^ 



dEhnepreenll 19921) . In addition to these two rates, there is 



third rate, the free-fall rate tjff ^ (Gpy^^, which is estimated 
from the local volume density. For an isothermal disk in hy- 
drostatic equilibrium, the midplane density can be derived 
from the scale height h = c^/iirCS) using po = S/(2/i), so it 
follows that ujff ^ u!j to within factors of a few. 

An approximate adjustment for the case of a two- 
component disk of gas and stars can be made by identify- 
ing the Je ans rate with the vertica l oscillation rate in the disk, 
shown bv lTalbot & Arnetti (Il975h to be: 



c* c. 
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where c* is the vertical stellar velocity dispersion. 



(2) 
(3) 



Similarly, the free-fall rate for gas at the equilibrium mid- 
plane density po can be written as 
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usi ng the appro ximation for the midpl ane gas den s ity de rived 
bv iLerov et all 12008.) from work by lElmegreenI (119891) and 
IKrumholz & McKed (5005l) : 
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where 77 = 1 + (cg/c*)(E*/Eg) is a correction factor that rep- 
resents the additional compression of the gas layer due to the 
gravity of the stellar disk. 

The importance of the correction factor 77 lies in its possible 
effect on the observed star formation law. For a star formation 
rate oc Yi^uji or oc SgO^ff, one expects Ssfr « for constant 
gas velocity dispersion c^. However, including the effect of 
the stellar disk increases the star formation rate by a factor 
of 77 (adopting the Jeans rate) or 77'/^ (adopting the free-fall 
rate). Since ?; generally anti-correlates with Sg, the tendency 
is to reduce the Schmidt power law below N=2, by an amount 
whic h dep ends on the relative dominance of the stellar disk 
(see 35.4l for further discussion). 

While the orbital rate is not directly tied to growth rates for 
gravitational instabilities, it may be more straightforward to 
estimate in numerical simulations since it reflects the grav- 
itational potential of the galaxy as a whole. In addition, if 
star formation feedback constrains the Q parameter towards a 
narrow range of values, whereas the assumption of constant 
Cg is not generally valid, then a star formation law based on 
the orbital rate may turn out to be more useful in practice, 
since 11/2 ^n^Y be easier to estimate than Sg/cg. On the other 
hand, if Q varies with radius, the orbital and Jeans rates may 
be substantially different. Note that when the stellar disk is 
taken into account, the ratio cjj.sg/ 17 oc = 2~' +2*' = Geff' 
which is approximately the instability parameter for a com- 
bined disk of gas and stars derived by Wang & Silk ( 1994) (ig- 
noring the anisotropy of the stellar velocity dispersion). Thus, 
a proportionality between ojj sg and 17 is maintained as long as 
2eff is roughly constant. 

Hereafter we omit the subscript 'sg' so that timescales are 
calculated including the correction factor for the stellar disk. 

3. OBSERVATIONAL DATA 

We primarily make use of radial profiles of E.^ , Esfr, and 
Eg [co mposed of E(H I) and E(H2)] published by Lerov et al] 
( I2OO8I hereafter L08), which we briefly describe here. The 
profiles of E^. are derive d from Spitzer 3.6 pm maps from 
SINGS (IKennicuttet all 12003) using an empirical scaling 
factor from K to 3.6 pm and a /iT-band mass-to-light ratio 
of Tf = Q.SMq/Lqk- Esfr is derived from a combina- 
tion of FUV images from the GALEX Nearby Galaxies Sur- 
vey dGil de Paz et alj|2007h and 24 ^m images from SINGS. 
E(H2) is derived from CO(2-l) images from the IRAM 
30m telescope assuming /co(2-l)//co(l-0)=0.8 and a Galac- 
tic value of the CO-to-Hi conversion factor, Xco=2 x 10^" 
cm-2 ( K kms"')"^ S(Hl) is derived from the VLA THINGS 
survey (lWalteretalJl2008[) assuming the 21 -cm line is op- 
tically thin. We choose three galaxies from L08's "dwarf" 
sample (NGC 925, 2403, and 2976) and six galaxies from 
their "spiral" sample (NGC 2841, 3521, 4736, 5055, 6946, 
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Fig. 1 . — Star formation an d dynamical timescal es plotted vs. radius normalized to Ris- (Left) Timescales are based on the radial gas and stellar profiles of LOS 
and the HI rotation curves of lde Blok et al.l 120081) . (Right) Same, except BIMA CO profiles are used for H2 and 2MASS A^5-band profiles for S,. Line colors 
indicate different timescales, line styles indicate different galaxies. 



and 7331). Radial profiles of the dwarf sample were obtained 
at resolutions of 9"-26" (chosen to achieve a spatial resolu- 
tion of 400 pc) while radial profiles of the spiral sample are at 
resolutions of 11 "-35" (chosen to achieve a spatial resolution 
of 800 pc). 

We u se the BIMA Sur vey of Nearby Galaxies (BIMA 
SONG; iHelferetalJ 12003 !) as an alternative source for CO 
data. The publicly available CO(l-O) data cubes for the 
nine sample galaxies were convolved to the resolution of the 
THINGS H I data (typically --10"), or, in cases where the H I 
data have higher resolution, both CO and H I were convolved 
to a resolution of 8". Radial profiles were derived by averag- 
ing in elliptical rings using the center positions, inclinations, 
and position angles of L08. H I profiles were derived from 
the integrated intensity maps provided on the THINGS web- 
site. CO profiles were derived from integrated intensity maps 
produced using a smooth-and-mask technique: a mask was 
generated using a 3a threshold in the cube after smoothing to 
15", and unmasked pixels exceeding la were summed in ve- 
locity. For galaxies with very weak CO emission (NGC 925, 
2403, 2841, and 2976), the mask was generated using a 6a 
threshold in the H I cube, again smoothed to 15". We verified 
that the profiles derived from the masked cubes were consis- 
tent (within the uncertainties) with profiles derived by direct 
summation of the channel maps. To derive the total (H I-1-H2) 
gas profile we extrapolated the CO profile beyond the last re- 
liably measured point (typically R=90") using an exponential 
fit to the region R > 40". 

To obtain a profile appropriate for the higher resolution 
of the SONG data, we use A"-band i mages from the 2MASS 
Large Galaxy Atlas d Jarrett et al.l2003f) . available on the IRSA 
website, assuming a native resolution of 2. "5 and applying a 
Gaussian smoothing kernel to achieve the same resolution as 
the CO and H I images. The radial profile is taken as the me- 
dian surface brightness in each elliptical annulus, corrected 



for inclination and converted to a mass surface density using a 
mass-to-light ratio of Tf = 0.5 Mq /Lq k for consistency with 
L08. For EsFR we simply adopt the profiles provided by L08. 
Note that these are at somewhat lower resolution than the CO, 
H I, and stellar profiles, and so comparisons of tsf with dy- 
namical timescales must be made with caution, particularly 
near the galaxy centers. 

In order to estimate dynamical timescales, we also require 
radial profiles fo r f^, Cg, and Cy. W e use the THINGS H I 
rotation curves of lde Blok et al.| (|2008) to derive i7, and adopt 
a constant value of 10 km s"' for Cg, similar to that adopted 
by L08. Lacking information about the radial variation of the 
stellar velocity dispersion, we follow standard practice (e.g., 
L08) by assuming a constant stellar scale height and an 
isothermal disk such that 

= i± . 

27rGS* 

Here we define /z* such that an isothermal disk has p(z) oc 
sech^(z /2/i,,), wh ich asymptotically approaches exp(-z//i*) at 
large z (Ivan der K ruit 1988). We estimate the scale height us- 
ing the empirical relation between radial and vertical scale- 
lengths of .Kregel et al. (2002), 




= 7.3±2.2, 



taking for the exponential K-hmd scalelengths given by the 
2MASS LGA (parameter a). We then use this assumed value 
of h^, along with the profiles to derive the radial variation in 
. While the assumption of a sech^ disk is questionable, with 
numerical calculations bv lBaneriee & Jo j (120071) suggesting 
that in the presence of a gaseous disk the vertical stellar profile 
is closer to exponential than sech^, adopting different vertical 
profiles has little effect on the relation between c* and for 
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Fig. 2. — Comparison of Jeans timescales (in red) calculated from both stellar and gas profiles (left) and gas profiles only {right) from LOS. 



TABLE 1 

Radial scalelengths for various timescales 

CO-faint sample CO-bright sample 
/(tsf)/R25 0.34 ±.08 0.48 ±.14 

/(tj)/R25 0.65 ± .16 0.40 ± .06 

'(Tff)/R25 0.79 ± .26 0.38 ± .09 

1(toa)/R25 0.76 ± .23 0.50 ± .13 



a given (Ivan der KruijfToSSh . 

4. RESULTS 
4.1. Radial variation in timescales 

Figure [Ha) shows the radial variation in the observed star 
formation timescale, tsf = 5]g/EsFR, in comparison with the 
three dynamical timescales that we consider, based on the 
LOS data. Timescales are derived from rates using the rela- 
tion r = Itt/uj, with Torb = 27r/il. FigurefUb) shows the same 
quantities but for the BIMA SONG data. Different line styles 
denote different galaxies, whereas different colors denote dif- 
ferent timescales. Radii have been normalized by the optical 
radius R25 as given by LEDA, typically w300". The CO-faint 
galaxies NGC 925, 2403, 2841, and 2976 are shown in the 
upper panels, while the CO-bright galaxies NGC 3521, 4736, 
5055, 6946, and 7331 are shown in the lower panels. Note 
that Tsf has been scaled by 0. 1 to make it comparable to the 
other timescales. 

Generally speaking, tsf shows an increase with radius, 
from ^1 Gyr in the central regions to ^10 Gyr at the edge 
of the optical disk. The central values of tsf appear to be 
smaller in CO-faint galaxies; as discussed by LOS, this pro- 
vides indirect evidence that the H2 content in these galaxies is 
underestimated, since it seems unlikely that CO-faint galax- 
ies would have higher star formation efficiencies, given their 
low metallicities, strong radiation fields, and weak gravita- 
tional fields. For the CO-bright galaxies, tsf appears to flat- 



ten in the inner regions; in these regions the ISM is largely 
molecular, and a near proportionality between Ssfr and is 
found (constant tsf). NGC 6946 is a notable departure from 
these overall trends: for this galaxy tsf reaches a minimum 
at intermediate radius of '^0.57?25, due to the very active star 
formation in its spiral arms. 

Table [T] summarizes the radial exponential scalelengths for 
the various timescales, normalized to R25. These are ob- 
tained from unweighted least-squares fits to the curves in Fig- 
ure The most notable trend is that, for the CO-faint 
galaxies, the exponential scalelength for tsf is a factor of 
^2 shorter than the scale lengths of all of the dynamical 
timescales. Thus none of the dynamical timescales seems able 
to predict the star formation timescale in these galaxies. For 
the CO-bright galaxies, on the other hand, the radial increase 
in Tsf is reasonably well matched by the radial increase in the 
dynamical timescales. 

Figures |2a) and|2b) show the effect of neglecting the stel- 
lar contribution to the Jeans rate when calculating tj. Partic- 
ularly in the CO-faint galaxies, where the radial gas profile is 
dominated by H I with fairly constant surface density, consid- 
ering only gas leads to a very flat profile for tj. In CO-bright 
galaxies, the gas-only tj profiles also show much more vari- 
ation with radius, as expected given the considerable amount 
of substructure in the radial gas profiles. 

Figures 13 a) and |3lb) display the same data as in Fig. [T] 
but now the ratio of tsf to each of the dynamical timescales 
is plotted. These plots reveal more clearly that in CO-faint 
galaxies, tsf displays a steeper radial gradient than do any of 
the dynamical timescales. For CO-bright galaxies, the ratios 
show little systematic variation with radius, except near small 
R where the flattening of tsf and the continued decline of tj 
and Tff is apparent. 

In general, the higher resolution SONG data [Figures [Tib) 
and lllb)] are consistent with the LOS profiles outside of the 
central regions (R/R25 > 0.1), as expected since the differ- 




Fig. 3. — Ratios of star formation to various dynamical timescales plotted vs. radius normalized to S25. Data are from LOS (left) and BIMA SONG (right) as 
in Fig.[T] Line colors indicate different timescales, line styles indicate different galaxies. 




Fig. 4. — Comparison of star formation timescale with (a) Jeans timescale for both stars and gas; (b) midplane free-fall timescale; (c) orbital timescale. Data 
are from LOS. The soHd line represents the average least-squares slope for the sample, while the dotted line represents a slope of unity. 



ence in resolution only becomes apparent on scales of <20". 
Closer to the center {R/R25 < 0.1), differences become ap- 
parent at high resolution. As discussed in ^5.11 these trends 
may reflect limitations in estimating gas masses or dy namical 
timescales in these regions. For the analysis in ^4.21 we con- 
fine our analysis to the LOS data, which are more complete in 
radial coverage and are also well-matched in resolution. 

4.2. Correlations Between Timescales 

Figure |4] shows correlation plots between tsf and each 
of the three dynamical timescales. These plots demonstrate 
clearly the minimum value of tsf (about 1-3 Gyr) achieved in 
the inner parts of the CO-bright galaxies (black, blue, and ma- 
genta curves). In the CO-faint galaxies, the slope of the rela- 
tion is generally flatter than unity, again reflecting the steeper 
radial gradient of tsf- Because of the scaling factor of 0.1 



applied to the abscissa, tsf is always longer than any of the 
dynamical timescales considered. 

Figure |5] shows similar correlation plots for the case in 
which only the H2 gas is used to calculate the star forma- 
tion, free-fall, and Jeans timescales. Such a comparison is 
motivated by the possibility that H I is inert as far as star for- 
mation is concerned, and that warm H I in particular may not 
contribute to the gravitational instability of the disk. We find 
no correlation between tsf.Hi and any of the three dynamical 
timescales — indeed, they appear to be anti-correlated, partic- 
ularly in NGC 4736 and 6946. Thus, if star formation only 
involves H2, there is no evidence that a dynamical timescale 
is relevant for establishing its rate. The negative slopes in the 
correlation plots, which at face value suggest that molecular 
clouds produce stars with greater efficiency in the outer disk, 
may indicate that CO emission provides an incomplete inven- 




Fig. 5. — Compaiison of star formation timescale with (a) Jeans timescale for both stars and gas; (b) midplane free-fall timescale; (c) orbital timescale. Data 
are from LOS, and only the H2 component of the gas is used to compute the timescales. The dotted line represents a slope of unity. 




Fig. 6. — Comparison of star formation timescale, using only the H I component for the gas mass, with (a) Jeans timescale for both stars and gas; (b) midplane 
free-fall timescale; (c) orbital timescale. Data are from LOS. The solid line represents the average least-squares slope for the sample, while the dotted line 
represents a slope of unity. 



tory of H2 mass in these regions. 

Finally, in Figure |6] only the H I component of Sg is used 
to calculate the star formation timescale, though the H2 is in- 
cluding in estimating the dynamical rates. If star formation in 
giant molecular clouds (GMCs) occurs at a constant rate per 
unit H2 mass, and if GMC lifetimes are not a function of en- 
vironment, then any variation in the star formation timescale 
(within or among galaxies) should reflect the timescale for 
forming GMCs from H I. Thus in a steady-state situation, 
the rate at which H I clouds form GMCs (e.g., the Jeans 
rate) should be proportional to the rate at which H I would 
be consumed by the current SFR. While all three dynamical 
timescales show a good correlation with tsf.hi, the correla- 
tions are much flatter than linear (power law slopes 0.4-0.5). 
In other words, regions of low SFE possess dynamical times 
that are too short, and regions of high SFE possess dynamical 
times that are too long, compared to the simple steady-state 
model. Of the three dynamical timescales, tj appears to show 
the strongest congelation with tsf.hi- 

5. DISCUSSION 
5.1. Behavior of Timescales Near Galaxy Centers 



The higher angular resolution of the BIMA SONG data 
[Fig-Eb) and Ob)] allows us to probe dynamical timescales 
closer to the centers of the galaxies, though we note that tsf 
is still calculated using the lower resolution Esfr profile. For 
several of the galaxies, such as NGC 4736, 5055, and 6946, 
both Tj and Tff show a steep decline for7?<30" {R/R25 < 0.1) 
in the high-resolution data, reflecting the steep increase in 
associated with the stellar bulge. While the bulge is an im- 
portant additional term in the gravitational potential, we have 
almost certainly overestimated its effect by including it in 
but adopting a disk-based value for c». A larger value of c» 
would serve to reduce uiff and ujj in these regions and thus 
increase the corresponding timescales. Moreover, the appli- 
cability of the Jeans time for a thin disk becomes increasingly 
suspect within the more spherical potential of a stellar bulge. 
Thus, the short dynamical timescales in this region may be 
more indicative of limitations in our modeling than of condi- 
tions that are especially conducive to star formation. 

The ratio of tsf to the dynamical timescales shown in 
Fig. |3lb) also shows the effect of the sharp central drop in 
Tj and Tff, but the ratio is enlarged even further by an appar- 
ent increase in tsf near the centers of several galaxies (again, 
NGC 4736, 5055, and 6946 are prime examples). Although 
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this may be due in part to the mismatched resolution of the 
SsFR profiles, these trends are also exhibited to a lesser ex- 
tent in the LOS data [Fig.lTfa)]. A large reservoir of molecular 
gas with relatively little massive star formation has been noted 
in the central region of NGC 4736 ( Wong & Blitz"2000), but 
other explanations exist as well, such as a change in the CO- 
to-H2 conversion factor Xco due to higher gas pressures (e.g., 
|Reganetal."2001; see also discussion in ^5.2l below). 

Notable in Fig.lUb) is the sharp drop in rsp/Tff in the central 
region of NGC 2841; this region exhibits very little neutral 
gas but a high rate of star formation, as seen in both FUV and 
24fim emission (LOS). The dearth of gas led to the omission 
of this region from the radial profiles of LOS, and hence it is 
also omitted from Figures [Tt a) and[3a). It is unclear whether 
this galaxy has experienced a recent starburst that has depleted 
the gas, or whether the central FUV and IR emission arises 
from a source other than star formation. 

5.2. Variations in the Ratio of Timescales 

The results of ^indicate that, ignoring any distinction be- 
tween H I and H2, the star formation timescale cannot be 
modeled as simply proportional to a dynamical timescale. In 
low-mass, CO-faint galaxies, the dynamical timescales fail to 
reproduce the steep radial gradient observed in rsp. In CO- 
bright galaxies they fare somewhat better, but are unable to 
reproduce the constant tsf found in the inner, H2 -dominated 
regions and also predict too shallow a radial gradient in the 
outer disk. These conclusions were akeady reached by LOS 
for Tff and Torb, so it is not surprising that they apply to tj as 
well. 

One possibility is that tsf is systematically overestimated, 
due to an overestimate of the amount of star-forming gas or an 
underestimate of the SFR. This would have to occur both in 
the inner, H2-dominated disks of spirals (where the discrep- 
ancy is 0.6-O.S dex) and in the outer disks of H I-dominated 
dwarfs (where the discrepancy is O.S-1 dex). In the inner 
disks, a decrease in Xco by a factor of 4-6, or an increase 
in the SFR by a similar factor, would help to linearize the re- 
lation between tj and rsp. There is some evidence for at least 
the first of these effects: a decrease in Xqo by factors of 2-5 
has been inferred towards galaxy centers based on measure- 
ments of '^CO emission (Pa glione et al . 2001). On the other 
hand, it seems unlikely that SFRs are substantially underes- 
timated by LOS, since they employed a combination of FUV 
and IR data. In outer disks, tsf could be overestimated if a 
large fraction (S4%) of the H I is inert as far as star formation 
is concerned (e.g., in the warm neutral phase), or if the stellar 
initial mass function (IMF) dis favors high-mass stars du e to 
a paucit y of high -mass clusters (IWeidner & Kroupall2005l but 
see also lElmegre en 2006). 

Alternatively, the dynamical timescale rj may be understi- 
mated, or equivalently cjj overestimated. In outer disk re- 
gions, this could be sympto matic of a dec line in the stellar 
mass-to-light ratio. While Bell et al.l (I2003h deduce a modest 
decrease \n M /Lk for bluer galaxy colors, the trend is highly 
dependent on metallicity, and is too small (<0.2 dex) to ac- 
count for the observed values of tj. 

Another source of systematic error may be our crude ap- 
proximations for the velocity dispersions Cg and c, . It seems 
unlikely that an increase in Cg or c» would occur in the outer 
disk (such trends have never been reported), but an unmod- 
eled increase of one or both of these might lead us to under- 
estimate Tj in inner, H2-dominated regions. As noted above 
(EB> c* may increase faster than expected near galaxy cen- 



ters due to the presence of a bulge. An increase in the H I 
velocity dispersio n towards the centers of THINGS galaxies 
has been found by'Tambu rro et al.l (120091) . though the increase 
appears too small (about a factor of 2) to offset the apparent 
decline in tj as /? ^ 0. 

In summary, while inclusion of a number of effects could 
act to linearize the relation between tsf and tj, the most plau- 
sible appear to be a change in Xqo and c, in the inner disk and 
a change in the fraction of star-forming gas in the outer disk. 
Investigating each of these possibilities remains an important, 
albeit difficult, observational challenge. 

5.3. Dual-Phase Approach 

The tighter correlation between tsf when only H I is in- 
cluded and Tj (Fig. |6]l suggests another possibility: that the 
dynamical timescale is important for the formation of H2, but 
that other physical processes (presumably internal to molecu- 
lar clouds) govern the timescale for star formation. The sim- 
plest model of this type would be one in which GMC lifetimes 
are constant and the SFR per unit H2 mass is also constant. In 
that case the star formation rate would be proportional to the 
molecular cloud formation rate (CFR): 



E(HI) 



Tj 



(6) 



However, rearranging this implies that tsf.hi « rj, which is 
not observed. 

A slightly less restrictive model keeps the SFR per unit H2 
mass constant, and assumes a steady-state cycling between 
HI and H2: 



I](HI) E(H2) 



■J SFR 



ThI^H2 TGMC,life TGMC,life 

Under the assumption that rni-tH, c>c tj, this implies that 



_ E(HI) 

TSF,HI = OC 



Es 



FR 



TGMClife 



(7) 



(8) 



The (approximately) observed relation rj cx (tsf.hi)"^ would 
then require 

TGMClife OC Tj"' . (9) 

Independent (since it does not rely on the SFR) support for 
the relation in Eq. ^ comes from the observed empirical cor- 
relation 

D — S(H2) 

Rmo\ = OC Ph , (10) 

where the midpl ane hydrostatic pressure /"h is given by 
(lElmegreenll989h : 



A: 



1 + 



2tt 



(11) 



The proportionality between Rmo{ and Ph, which holds approx- 
imately in the CO-bright galaxies exam ined by Wong & Blit3 
(120021) and lBlitz & RosolowskyI (12004 . implies that 



„ T'GMC.life 1 
"mol = OC 



(12) 



assuming that Cg and ?/ do not vary significantly with radius. 
Thus, if thi^h2 OC tj, then TGMCiife oc rf' . 

Indeed, the tightness of the relation between R^o\ and Pj, is 
directly related to the tightness of the relations between rj and 
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Fig. 7. — Time evolution of Eg, rgp, and tlie instantaneous Schmidt law 
index N = dInSspR/rflnEj for a closed-box model (solid line) where the 
star formation timescale is proportional to the Jeans time, and Cg/c, is fixed 
at 0.2. The gas surface density and star formation timescale are normalized 
to their initial values, and tq is the initial star formation timescale. Constant 
Schmidt law indices of 1 (magenta dotted line), 1.5 (blue short-dashed line), 
and 2 (red long-dashed line) are shown for comparison. 

TsF,Hi and between Ssfr and S(H2). If Ssfr « S(H2) then it 
follows that TsF Hi oc R^„\, at the same time, the definition of 
Ph implies that rj oc P^^, setting aside the correction factor rj. 



This may provide a route for understanding the origin of the 
^moi oc Ph relation in terms of molecular cloud formation and 
destruction, and for unders tanding the limitation s to its wider 
applicability. For instance, 'Koyama & Ostrike^ ( l2009b ) have 
found that in simulations where the Q parameter is allowed to 
vary, the /?moi « Pj, relation no longer holds. This suggests that 
galactic rotation has a strong influence on thi^Hi or TcMCiife 
which is obscured by the fact that ojj w 17 in real galaxies. 

Thus, one can account for the two empirical correlations, 
SsFR OC S(H2) and tj oc (tsf.hi)"'^, by assuming that H I forms 
H2 on a Jeans timescale, and that the reverse process of disso- 
ciating H2 occurs on a timescale that scales inversely with the 
Jeans timescale. In principle this could be a viable scenario, 
since conditions that favor rapid molecular cloud formation 
(such as high external pressure, or a high gas shielding col- 
umn) could also forestall cloud destruction. For instance, if 
the expansion of H II regions is related to the destruction of 
molecular clouds, then since the stall radius scales inversely 
with external pressure P, cloud lifetimes may therefore scale 
with P. However, such a scaling would be difficult to rec- 
oncile with the expectation that high-mass clouds are more 
quickly unboun d by the massive stars which form within them 
jMatzn er 2002). Further work will be needed to cast Eq. (|9]l 
in a dimensionally correct form and to derive the observed 
proportionality constant. 

5.4. Comparison with the Schmidt Law 

Although we have not established a linear relationship be- 
tween the observed star formation timescale and the modi- 
fied Jeans timescale tj, it is nonetheless interesting to con- 
sider the implications such a relationship would have on the 
observed star formation law. For simplicity we consider a 



closed-box model, where all the gas is in place at f=0. By 
assumption the total mass surface density at a given radius is 
constant: S(r,r) = So('")- Defining the gas fraction as fg(r,t) = 
I]j(r,f)/S()(r), simple integration of the classical Schmidt law 
given by Eq. ([T]i yields 



f^=l [l+ktiN-DiEof-^]' 
] exp(-fef) 



1/(^-1) 



N=l 



fcf. lBell & de Jon^l2000h . In terms of the local star formation 
timescale, 

So 1 
'^0= — 



^sfr.o 



this can be written as: 



[l + (A^-l)f/To]" 
exp(-r/To) 



1/(^-1) 



N=l 



(13) 



If instead we adopt a star formation law based on the mod- 
ified Jeans rate oj], such that 



^SFR OC I iH — 

Cg \ C, l^g 



then once again starting from a pure gas disk, 



dt 



r 

To 



1+- 



(14) 



(15) 



IS quite 
1.5 Schmidt law. 



We have integrated this equation numerically for to=1 Gyr 
and a fixed ratio Cg/c*=0.2. Figure|7]shows the time evolution 
of the gas and star formation surface densities, normalized 
to their initial values. Note that the evolution of E 
similar to what would be expected for an A^: 
However, the measured Schmidt law index varies with time, 
from an initial value close to 2 (appropriate for a pure gas 
disk) to a value close to 1 at late times (when the gas is almost 
entirely consumed). This simple example illustrates that the 
observed Schmidt law index may be a function of the current 
evolutionary state of a galaxy. 

5.5. Relevance of Timescales Estimated from Average 
Densities 

The average gas densities used to derive tj and ts in this pa- 
per are often much less than the surface and volume densities 
of giant molecular clouds (^^100 Mq pc"^ and ^^100 cm"^ re- 
spectively) where the bulk of star formation actually occurs. 
One may certainly question whether such average densities 
can be reliably used to infer the SFR. We adopt this as a work- 
ing hypothesis, assuming that either (1) large-scale average 
densities correlate well with local GMC densities, as would 
be the case if both scale with ambient pressure; or (2) the star 
formation rate is ultimately limited by the GMC formation 
rate, which reflects the conditions of the ISM on much larger 
scales. The fact that Qeff ~ 1 in observed disks offers indirect 
evidence that large-scale dynamics influences star formation, 
since it is presumably star formation feedback that increases 
Cg and thus Q to values that are marginally unstable. Unfor- 
tunately, the GMC formation rate is difficult to infer directly, 
since the ages of clouds cannot be easily determined, and the 
overall HI-H2 balance depends on both H2 formation and de- 
struction rates. Detailed correlation studies of young stellar 
objects and molecular clouds in the nearest galaxies, such as 
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the Large Magellanic Cloud, should be able to better constrain 
these parameters. 

Nonetheless, we cannot exclude the possibility that the 
star formation and dynamical timescales may not scale lin- 
early because they are determ ined over regions w i th very 
different densities. Recently, iKo yama & Ostrike rl (I2009al) 
have conducted numerical simulations of galactic disks in- 
cluding stellar feedback, and compared the large-scale dy- 
namical timescales (Tff, Torb, and rj gas) with star formation 
timescales derived by assuming a constant star formation ef- 
ficiency above a density threshold. They also find a non- 
linearity between tsf and Tff, which they attribute to estimat- 
ing the dynamical timescales on spatial scales much larger 
than star-forming clouds. In particular, tsf / Tff increases in re- 
gions of high Ey, similar to what we observe in H2-dominated 
regions. Intuitively this makes sense, since tsf becomes con- 
stant as the threshold density is approached, whereas Tff con- 
tinues to decrease. On the other hand, it is not clear whether 
such simulations can account for the minimum in the tsf /Tff 
ratio at intermediate radii and its subsequent increase in the 
outer, H I-dominated disk. 

Finally, we note that any theory relating dynamical 
timescales with tsf must explain the order of magnitude dif- 
ference between them. The difficulty in accounting for this 
low efficienc y factor has been d iscussed extensively in the lit- 
erature (e.g., lElmegreenll2002l) . One possibility is that this 
factor arises from a universal density structure generated by 
turbulence, coupled with a threshold density for star forma- 
tion that is a fixed multiple of the mean density. This ef- 
fectively ensures that once H2 forms, only a small, roughly 
constant fraction of it participates in star formation. Varia- 
tio ns of this idea have b een explored by 'Elmeg reenI (l2002h 
and lKrumholz & McKeel (2005), but the appropriate value for 
the threshold density has yet to be clearly established. 

6. CONCLUSIONS 

We have determined the dynamical timescales in nearby 
disk galaxies using three common estimators, modified to ac- 
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count for the gravity of the stellar disk. The Jeans and orbital 
timescales are related by an effective Q parameter, which is 
observed to be close to 1 in real galaxies, making the two 
timescales interchangeable. Use of the modified Jeans rate to 
derive the SFR leads to a star formation law that resembles 
the A^=1.5 Schmidt law, but which, for a closed-box model, 
has a value of which declines over time. On the other hand, 
recent observational data suggest that the star formation rate 
is not simply proportional to the modified Jeans rate, although 
inclusion of the stellar component does improve the correla- 
tion. Several effects, including a change in the Xqo factor 
in the inner disk, improper modeling of the stellar velocity 
dispersion of the bulge, and a decrease in the fraction of star- 
forming gas in the outer disk, could contribute to the non- 
linearity. Future studies using the far-infrared dust emission 
to derive Xco and resolving the cloudy structure of the neutral 
ISM may be able to quantify these effects. 

Alternatively, we consider a two-step model where the 
Jeans rate determines the GMC formation rate, whereas 
GMCs form stars at a constant rate. In agreement with this 
model, we find a much tighter correlation between the H I 
depletion timescale (due to star formation) and the Jeans 
timescale, but again it is not linear; to linearize it would re- 
quire that the GMC (or H2) destruction rate is inversely cor- 
related with the Jeans rate. Although somewhat unexpected, 
such a scaling would also account for the correlation between 
the H2/H I ratio and the midplane hydrostatic pressure. Ob- 
servations of GMCs and their embedded stellar populations in 
Local Group galaxies will be needed to investigate this possi- 
bility. 
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